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Can superhorizon perturbations drive the acceleration of the Universe? 
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It has recently been suggested that the acceleration of the Universe can be explained as the 
backreaction effect of superhorizon perturbations using second order perturbation theory. If this 
mechanism is correct, it should also apply to a hypothetical, gedanken universe in which the subhori¬ 
zon perturbations are absent. In such a gedanken universe it is possible to compute the deceleration 
parameter qo measured by comoving observers using local covariant Taylor expansions rather than 
using second order perturbation theory. The result indicates that second order corrections to qo 
are present, but shows that if qo is negative then its magnitude is constrained to be less than or of 
the order of the square of the peculiar velocity on Hubble scales today. We argue that since this 
quantity is constrained by observations to be small compared to unity, superhorizon perturbations 
cannot be responsible for the acceleration of the Universe. 

PACS numbers: 98.80.Cq 


The observed acceleration of the expansion of the Uni¬ 
verse Qi is a profound mystery. It is usually explained 
by positing a new form of matter with negative pressure, 
so-called dark energy [sj, or by a modification of gen¬ 
eral relativity at large distance scales (4]. Recently it has 
been suggested that the acceleration is instead driven by 
the backreaction effect of inflation-generated superhori¬ 
zon perturbations via second order perturbation theory 
[^,|g . The backreaction of perturbations had earlier been 
studied in other contexts by Brandenberger and collabo¬ 
rators Q and others [g]. 

The basic idea is the following. As is well known, scalar 
perturbation modes whose wavelengths today are smaller 
than A ~ 8 Mpc have entered the nonlinear regime and 
are responsible for galaxies and galaxy clusters, while 
longer wavelength modes and in particular superhorizon 
modes (A > 3 Gpc) are still in the linear regime today. 
That is, the fractional density perturbation due to these 
modes is small compared to unity, so each individual 
mode evolves with high accuracy according to the lin¬ 
earized equations. Nevertheless the net effect of the back- 
reaction from all the superhorizon modes can still be sig¬ 
nificant. If we denote by e = 6p/p ~ 10 -4 the fractional 
density perturbation on Hubble scales today, then one 
naively expects second order corrections to be of order 
e 2 ~ 1CU 8 which is negligible. However, a more refined 
estimate for some particular second order effects obtained 
from integrating over all the superhorizon modes Q gives 
the scaling e 2 F{k la i n , fc max ), where fc m ; n and fc max are the 
minimum and maximum comoving wavenumbers of the 
scalar perturbation spectrum integrated over. Here F 
is a function for which F —> oo as fc m ; n —* 0. Thus, 
if the spectrum extends to sufficiently large wavelengths 
(as would be generated by a sufficiently long period of in¬ 
flation), second order effects can become significant. In 
particular, Refs. [ 1 , 6 ] argue that such effects can give 
rise to an effective cosmological constant and drive the 
present acceleration of the Universe. 


A potential difficulty with this idea is that it appears to 
be in conflict with the locality and causality of general rel¬ 
ativity. Specifically, consider the finite spacelike hyper¬ 
surface V given by the interior of the intersection of our 
past lightcone with the spacelike hypersurface of some 
fixed redshift z = zq. We will take Zo = 5 say so that 
the observed supernova used as evidence of the universe’s 
acceleration are to the future of V. Then, the initial data 
within V are sufficient to determine all the observations 
we make; in particular the observable effect of superhori¬ 
zon perturbations must be encoded in this initial data 1 . 
In addition observations strongly constrain this initial 
data: we know that the geometry of V on large scales can 
be accurately modeled as a Friedman-Robertson-Walker 
(FRW) background plus fractional perturbations of or¬ 
der e ~ 10 -4 . Thus the issue is whether or not the space 
of initial data for the gravitational and matter fields on 
V contains enough freedom, given the observational con¬ 
straints, to mimic the effects of dark energy when no such 
dark energy is present in the matter stress energy tensor. 

In this paper we shall argue that the freedom is in¬ 
sufficient. More precisely, we focus on the deceleration 
parameter q 0 = — a(t)a(t) / a(t) 2 , where a{t) is the scale 
factor, whose measured value is qo ~ —0.5. We will ar¬ 
gue that (i) The local spatial curvature within V is uncon¬ 
strained and can be altered by second order effects, giving 
rise to changes in qo that could in principle be of order 
unity. However this effect cannot change a positive value 
of qo to a negative value of qo- (ii) Non-isotropy and non¬ 
homogeneity of the initial data in V on large scales can 
give rise to negative values of qo, as suggested by Refs. @. 


1 Note that this initial data does not itself contain superhorizon 
modes: V is a finite spherical region whose comoving radius is ~ 
1.8 Hq 1 for zq = 5. So for modes which are super horizon today, 
V extends across only a small fraction of the mode wavelength 
at z = zq. 
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However the magnitude of this effect is constrained to be 
of order the square of the velocity perturbation on Hub¬ 
ble scales, which is constrained by observations of low 
order multipoles of the cosmic microwave background to 
be small compared to unity. 

Method of analysis: If the superhorizon perturbation 
mechanism for driving acceleration is correct, it should 
apply not just to our Universe but also to other hypo¬ 
thetical universes. For ease of analysis, we will analyze a 
fictitious, gedanken universe which differs from ours only 
in that the perturbation spectrum at early times is mod¬ 
ified to suppress the perturbations which are subhorizon 
today. In this universe the subhorizon perturbations are 
negligible today, while the superhorizon perturbations 
are taken to be the same as those used in Refs, js, 0- A 
demonstration that the superhorizon perturbation mech¬ 
anism does not work in this context is fairly strong ev¬ 
idence that it cannot work in our Universe. The only 
possible loophole is the possibility that subhorizon per¬ 
turbations somehow play an important role, which does 
not seem to be indicated by the analyses of Refs. [1, :6] 2 . 

In this gedanken universe, the length and time scales 
over which the gravitational and matter fields are varying 
are all of order Hq 1 or larger. This allows us to perform 
an analysis in a local region using Taylor series expan¬ 
sions of the Einstein and hydrodynamic equations, which 
is much simpler than second order perturbation theory 
about a FRW background. 

We model the matter source by the fluid stress en¬ 
ergy tensor T a0 = (p + p)u a u 0 + pg a0 , where p , p, u a 
and g a p are the density, pressure, 4-velocity and metric. 
Consider a comoving observer at an event V. Such an 
observer can measure the redshift 2 and luminosity dis¬ 
tance C of nearby events, and thus measure the redshift- 
luminosity distance relation C = C(z,9,<p). Here 9 and 
ip are spherical polar coordinates in the observer’s local 
Lorentz frame. The dependence on these angles arises 
since we are allowing general local solutions of the Ein¬ 
stein equations; there is no requirement of isotropy. For 
small 2 this relation can be expanded as Q 

C = A(9,p)z + B(9,f)z 2 + 0{z 3 ). (1) 

We define the Hubble constant H 0 and deceleration pa¬ 
rameter <70 as measured by the observer by comparison 
with the conventional FRW relation C = Hq 1 z+Hq 1 (1 — 
qo)z 2 /2 + 0(z 3 ), as in Ref. jfj: 

H 0 = (A - 1 ), q 0 = l~ 2 Hq 2 (A~ 3 B). (2) 

Here the angular brackets denote an average over the 
angles 0 , (p. 


2 In addition momentum conservation k « k\ + &2 rules out second 

order corrections to low spatial frequency observables (k < Hq) 
from interactions between very subhorizon modes k\ Ho and 
superhorizon modes &2 < Hq. 


The particular prescription m for angular averaging 
is chosen for later convenience. Note that there is no 
unique prescription; one could for example use the defi¬ 
nitions Hq 1 = (H) and q 0 = 1 — 2 H 0 (B). Observations 
that measure Hq and qo typically assume isotropy and 
therefore effectively angle-average at some stage of the 
analysis, but the precise nature of the averaging is not 
usually discussed. We will argue below however that the 
effect of this ambiguity on the values of H 0 and go is 
small. 

Using local Taylor series expansions we can compute 
Hq and go in terms of the the density, 4-velocity and 
velocity gradients of the cosmological fluid evaluated at 
the observer’s location V. We decompose the gradient in 
the usual way as 

^a^/3 — 7^(l?a/3 T u a up) T C a0 T io a0 u a ap, (3) 

where 9 , cr a0 , uj a0 and a a are the expansion, shear, vor- 
ticity and 4-acceleration. For Hq we find the well-known 
result 

Ho = (4) 

the locally measured Hubble constant is just the expan¬ 
sion of the fluid. For q 0 we obtain 


47r 1 

® = 3^ ( '’ + 3p) + 3JfS 

a(3 
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Now to a good approximation in the present epoch we 
have p = 0 (assuming cold dark matter and baryons with 
no dark energy), which implies from V a T a P = 0 that 
a a = 0. This yields 


Qo 


47T | 1 

3H% P+ 3Hl 


7 


a/3 





( 6 ) 


Discussion: Consider first the first term in Eq. ©. This 
reduces to the conventional value go = 1/2 for a flat, 
matter dominated Universe when Hq = 8np/3. However 
in the present context the relation Hq = 8np/3 need not 
be satisfied; H 0 is instead given by Eq. 0 . The deviation 
of this first term from 1/2 is related to the fact that the 
local analysis allows spatial curvature. If we define an 
effective local flk by f \ = 1 — 8tt p/ (3Hq) , then we obtain 
for the first term go = (1—Hfc)/2, the conventional answer 
for a Universe with matter and spatial curvature. 

The key point about the first term in is that it is 
positive. Hence this term cannot drive an acceleration. 

Consider next the second and third terms in Eq. ®, 
the squared shear and squared vorticity. These quan¬ 
tities have an unambiguous operational meaning; they 
can be measured by the observer in her local Lorentz 
frame. We can estimate the sizes of these terms as 
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&ab& ab ,u a bu ab ~ {Sv) 2 /l 2 , where Sv is the typical scale of 
peculiar velocity (deviation from uniform Hubble flow), 
and l is the lengthscale over which the velocity varies. 
In the present context we have l > H 0 , by our assump¬ 
tion that subhorizon modes are negligible, which implies 
that the contribution from the second and third terms in 
© to qo is of order Sqo ~ {Sv) 2 ~ £ ~ 10 4 . We con¬ 
clude that it is impossible in this model to achieve the 
measured value qo ~ —0.5 of the deceleration parameter. 

Note that the key difference between our analysis and 
that of Refs. Hi is one of interpretation. Refs, 
predict changes to qo that are quadratic in the the first 
order perturbation variables, in agreement with our Eq. 
©. The new information provided by our analysis is 
that the quadratic terms are in fact locally measurable 
and represent degrees of freedom of the cosmological fluid 
rather than of the gravitational field. In the argument 
above, Sv characterizes the total deviation of the fluid 
velocity from an FRW background, including both first 
and second order perturbations. The contribution of the 
quadratic terms in Eq. © to qo are constrained to be 
small since observations constrain the total velocity per¬ 
turbation Sv. Thus, while an order-unity renormalization 
of qo from second order effects is possible in principle, our 
analysis implies that such a renormalization would also 
require second order contributions to the fluid velocity 
that violate observational bounds. 

Details of derivation: We use the local covariant expan¬ 
sion formalism based on bitensors E3- We denote by 
x a the coordinates of the event V where the observer is 
making observations, and by x a the coordinates of an 
event Q in the observer’s vicinity. We shall mostly be in¬ 
terested in the case where Q is on the past lightcone of V. 
We denote by A an affine parameter along the geodesic 
x a = z a { A) that joins Q and P, chosen so that A = 0 
at Q and A = 1 at V. We define Synge’s world function 
(the squared geodesic interval) via 

1 f 1 d7 a dz& 

a(x,x') = -j o dX g a/ 3 [z a (A)] (A) (A). (7) 


be expanded in a local covariant Taylor series as 

u a {x,x') = u a {x) + u af3 {x)a- a {x, x')a-p{x, x') 

+—u al3l {x)a- a {x, x')a-f 3 {x , x')cr ;l (x, x') 

+0{s A ). (9) 

Using Eq. (B14) of Ref. 0 we can evaluate the coef¬ 
ficients to give u a/3 {x) = — V^“u^(ar) and u a/3ry {x) = 
V(“V^tU)( x). 

The measured redshift of the event Q is given by the 
ratio of the inner products k ■ u evaluated at Q and V. 
Using the definition of u a this can be written as 


1 + z = 


( 10 ) 


Using the expansion © and the definition © this can 
be written as 

Z = U a,3 k a kf}S — t a/37 fc ce fc ( gfc 7 S 2 + 0(s 3 ). (11) 


Next we evaluate the luminosity distance C. This is de¬ 
fined so that 47 t£ 2 is the ratio between an energy emitted 
per unit time isotropically at Q and an energy per unit 
time per unit proper area received at V: 




( 12 ) 


These quantities can be evaluated using the geometric 
optics approximation to the scalar (or Maxwell) wave 
equation 6]. The stress-tensor for the radiation held is 
T a0 = A 2 l a lp , where l a is defined as being the set of 
null vectors at Q normalized according to l a 'U a = — 1, 
and then extended along Q’s future light cone using the 
geodesic equation. We define the affine parameter A by 
l = d/dX. The normalization conditions for the vectors 
k and l imply that l = k/{ 1 + z), and hence the affine 
parameters s and A are related by s = A/(l + z). The 
amplitude A satisfies the differential equation 


d{\nA)/dX = —0/2, (13) 


Then er ;Q ,( x,x') = V a <r(x,x') is the tangent to the 
geodesic at V. We define s(x,x') = —a ;a (x,x')u a (x); 
we will use s as our expansion parameter. We define the 
vector k a by 

a- a {x,x') = s{x,x')k a {x); (8) 

this is a future directed tangent to the geodesic which is 
normalized so that k a u a = — 1 at V. We define gff (x, x') 
to be the linear operator of parallel transport from the 
tangent space at Q to the tangent space at V, and we 
define u a (x,x') = g aa {x,x')u a ’{x'). This quantity can 


where 0 = V Q Z“ is the expansion. We choose the nor¬ 
malization of A so that H«l/AforA—>0 near Q. 

The energy flux at V can now be computed as 
dE/{dtd 2 A ) = T a pu a uP = A 2 {k a u a ) 2 = A 2 { 1 + z)~ 2 . 
The luminosity at Q can be evaluated by integrating 
the energy flux over a small sphere about Q of radius 
A; this gives dE/dt = A 2 {l a ru a ) 2 (47tA 2 ) = 47r. Combin¬ 
ing these results yields C = (1 + z)/A. Using the relation 
s = A/(l + z) we can rewrite this as 

£= {l + z) 2 sA{x,x')~ 1/2 , (14) 

where we have defined A = H 2 A 2 . This quantity satis¬ 
fies A —> 1 as Q —> V and also from Eq. m satisfies 
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the differential equation din A/(dX) = 2/A — 0. By com¬ 
paring with Eq. (32) of Ref. 42.1 we see that A(x,x') is 
the van Vleck determinant[lCl I12l|. Using the expansion 
A(x,x') = 1 + 0(s 2 ) given in Ref. E3 and combining 
Eqs. lllll and 11411 gives a relation between redshift 2 and 
luminosity distance C of the form 09 > where the coeffi¬ 
cients are 


A(d,tp) 

B(9,<p) 


(X7 a uP)k a kfi ’ 

2 (V“VV)fc a fc^ 7 

(\7 a U,P)k a kp + 2[(V“Ztd)/c a /c ( g] 3 ' 1 ’ 


We now evaluate the averages over angles. Using the 
definition Hq = (A -1 ) together with (k a ky 3 ) = (gap + 
Au a up)/3 yields the result ©, using (V a u l3 )uaUp = 
a a u a = 0. Note that if we use the alternative angle¬ 
averaging definition Hq = (A) we instead obtain H// 1 = 
((0/3 + aijUirij + din,) -1 ), where k a = (l,n l ) in the lo¬ 
cal comoving frame at V. Evaluating this average treat¬ 
ing the shear and acceleration as small compared to the 
expansion yields H 0 = 0/3 — 2a^aij/(50) — a^i/O + 
O(a 4 /0 3 ) + 0[(aiCLi) 2 /O 3 ]. Thus the different averaging 
prescriptions give different answers. However the frac¬ 
tional differences are of order a 2 /0 2 , which we have ar¬ 
gued above is of order e and is small. 

We now evaluate the angular average of the quantity 
A~ 3 B. We write this as (A~ 3 B) = J/2 + J, where 
/ = ((S7 a W^u^)k a kpky) and J = 2((V“/fc a ^) 2 ). For 
J we obtain J = 2((0/3 + aijUiTij + aiUj) 2 ) = 20 2 /9 + 
4<7y(7y/15 + 2oiOi/3. Using the formula (k a kpkj) = 
(gapUry + ga-yU /3 + g^ua)/3 + 2UaUpu-y, we obtain for 
/ the formula 


31 = u 7 V 7 V Q u“ + u 7 V Q V 7 u a + w q V 7 V 7 m“ 

+6 u a upu 1 S7 a S/ 0 u 1 . (17) 

We can rewrite the first term by commuting the covari¬ 
ant derivatives, which gives zz 7 Vq,V 7 m“ — R a pu a u 0 = 
V Q a“ — (V ct zi 7 )V 7 n“ — R a pu a u 0 . Using the decompo¬ 
sition this term can be written as V 7 a 7 — 0 2 /3 — 
<x a po' a0 + Lo a pu a0 — R a pu a u 0 . Similarly the second term 
in Eq. evaluates to V 7 a 7 — 0 2 /3 — cr a pa a0 +to a puj a0 . 
By differentiating twice the identity u a u a = — 1, the 
third term can be written as — {S7 a up)S7 a u 0 , which us¬ 
ing the decomposition 0 evaluates to —6 2 /3 — a a pa a0 — 
u a p(jJ a0 + a a a a . Finally a similar manipulation of the 
fourth term shows that it reduces to —6 a a a a . Combining 
these results and using the Einstein equation to replace 
Rapu a u^ with 47r(p + 3 p) gives 

I = ^V Q a“ - i 0 2 - a a pa a0 + \x> a pu a0 

~a a a a - y(p + 3p). (18) 

Now combining the results for / and J and substituting 
into the formula © for qo gives the result ©. 


Finally we note that using the alternative angular av¬ 
eraging definition q$ = 1 — 2Hq(B) would not change our 
conclusions. This average can be evaluated by expanding 
the denominators in Eq. m treating the last three terms 
in Eq. © as small compared to the expansion term, 
using the identity (. k a kpk^ksk s ) = 16zz a 'U3'u 7 ti ( 5M e /3 + 
16fl , ( a /3U 7 u (5 u £ )/3 + 5(a/35 7 5W), and performing manip¬ 
ulations similar to those used above. The modifica¬ 
tions to Eq. 0 that result are: (i) changes to the nu¬ 
merical coefficients of the shear squared and vorticity 
squared terms; (ii) the addition of a term proportional 
to H// 3 u a X7 a (crprya 0ri ) which is of the same order as the 
shear squared term; and (iii) the addition of terms that 
are smaller than the terms retained by one or more pow¬ 
ers of the small parameters a a pcr a0 /0 or \/u> a puj a P/0. 
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